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Abstract. 

We present analytic approximations of chiral SU (3) amplitudes for the extrapolation of lattice 
data to the physical masses and the determination of Next-to-Next-to-Leading-Order low-energy 
constants. Lattice data for the ratio Fk/F x is used to test the approximation proposed 1 . 
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INTRODUCTION 

Quantum Chromodynamics (QCD) is believed to be the correct Quantum Field Theory 
for the strong interactions in the Standard Model with quark and gluons as degrees 
of freedom. However, at low energies, due to the strong coupling of QCD, quarks 
do not show up individually but in bound states called hadrons. In the regime where 
the light hadron spectrum is observed, a perturbative approach is not feasible. Many 
efforts have been done to study this sector in the past years. Chiral Perturbation Theory 
[2, 3] as the correct Effective Field Theories of the Standard Model at low energies 
and simulations of QCD performed numerically on a finite lattice discretizing space and 
time are the most used ones. In particular, lattice simulations have made huge progress 
in the light quark sector employing quark masses corresponding to pion masses as low 
as 150 — 200MeV or even physical masses (see, for example, [4]). In case the simulation 
does not use physical values for the masses one needs a certain extrapolation to the 
physical point, extrapolation that is performed nowadays in different ways. The most 
elaborated once relies on Chiral Perturbation Theory (ChPT) which provides the correct 
analytic structure of amplitudes in terms of several a priori undetermined constants, 
the so-called low-energy constants (LECs), independent of the light quark masses by 
construction. The number of LECs increases with the order of the expansion in ChPT 
also increases. It is therefore impossible to extract all of them from experimental data. 
Lattice calculations offer a new scenario for determining LECs since one can tune the 
quark masses in a lattice simulation. State-of-the-art lattice studies use next-to-leading- 
order (NLO) ChPT for chiral extrapolations, determining for example several LECs at 
this & (p 4 ) order. NNLO ChPT results have only been used recently for the interpretation 



1 This work is based on Ref.[l] 
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FIGURE 1. Skeleton diagrams for the generating functional Z(, of ff(P 6 ). Simple dots, crossed circles, 
black box denote vertices from leading-order, NLO, NNLO Lagrangians respectively. Propagators and 
vertices carry the full tree structure associated with the lowest-order Lagrangian. 



of lattice data [5,6, 7]. 

For comparison with lattice calculations, the explicit dependence on the masses 
should be known. We propose analytic approximations of chiral SU (3) amplitudes for 
the extrapolation of lattice data to the physical masses and to determine NLO and NNLO 
LECs through a more interesting approach than a polynomial. 



A DISSECTION OF CHIRAL PERTURBATION THEORY AT 0{p 6 ) 

The most compact representation of ChPT in the meson sector is in terms of the gen- 
erating functional of Green functions Z [2, 3]. Analogous to the chiral Lagrangian, the 
generating functional permits a systematic chiral expansion: 

Z = Z 2 + Z 4 +Z 6 + ... (1) 

The NNLO functional Z 6 of 0(p 6 ) is itself a sum of different contributions shown 
pictorially in Fig. 1. In addition to tree diagrams of 0(p 6 ) (diagram g), there are two 
classes of contributions requiring separate treatments: irreducible (diagrams a,b,d) and 
reducible (diagrams c,e,f) contributions. For simplicity, we expose here the irreducible 
case (the reducible one can be found in [1]). With dimensional regularization, the irre- 
ducible diagrams have both double- and single-pole divergences. Moreover, the single- 
pole divergences of each irreducible diagram are in general non-local. Renormalization 
theory guarantees, however, that the sum of the three diagrams has only local diver- 
gences [8, 9], i.e., polynomials in momenta and masses in momentum space. Chiral 
symmetry guarantees that these divergences can be absorbed by the LECs of 0(p ) via 
diagram g. In this process an arbitrary renormalization scale /i is generated. The sum of 
diagrams a,b,d,g is then finite and can be written in the form 
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The monomials O a (x) (a = 1, . . . ,94) define the chiral Lagrangian of 0(p 6 ) [10] with 
associated renormalized LECs C£(ju), the (i = 1, . . . , 10) are renormalized LECs 

of 0(p 4 ) with associated beta functions T ; - [3] and the coefficients and rj^ are 

listed in Ref. [9]. Fq is the meson decay constant in the chiral SU (3) limit. The chiral log 
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involves an additional (arbitrary) scale M but Zg~ ru ~ ru g as we n as th e total generating 
functional are independent of both \x and M. Hi(x;M) are one-loop functionals 
associated with diagram d whereas the two-loop contributions (except for the chiral 
logs) are contained in the functional K(x;M). The functional (2) is scale independent 
[1]. 

The complete generating functional of 0(p 6 ) is then given by the sum of the irre- 
ducible and reducible parts, again independent of both scales /i and M: 
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Analytic Approximation for chiral ST/ (3) 

The genuine two-loop contributions contained in the functional K(x;M) are usually 
only available in numerical form for chiral SU(3). On the other hand, the one-loop 
contributions can be given in analytic form and the dependence on meson masses 
is manifest. For the chiral extrapolation of lattice results, we therefore suggest the 
following approximate form for the functional of 0(p 6 ): 
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The approximation proposed, called Approximation I from now on, consists on drop- 
ping the K(x;M) pieces and the reducible 2-loop contributions from diagrams a and c in 
Fig.l while keeping all the chiral logarithms L (more details can be found in Ref. [1]). 
For a reliable determination of renormalized LECs the analytic approximation (5) should 
be scale independent as it is the case. The main issue to address now is the applicability 



TABLE 1. Summary of results for the fitted Fk/F k for both Approximation I 
and Approximation II. 





L 5 • 10 3 


(q 4 +q 5 )-io 3 Gev 2 


(q 5 +2q 7 )-io 3 GeV 2 


Approx I 
Approx II 


0.76 ±0.09 
0.75 ±0.09 


0.37 ±0.08 
0.20 ±0.07 


1.29±0.16 
0.71 ±0.15 



of that approximation. The simple answer is based in a naive chiral counting. Consider- 
ing a SU (3) amplitude normalized to 1 at LO, the chiral counting suggest that the NLO 
to be of the order of 0.3, the NNLO of the order 0.3 2 = 0.09 and the NNNLO of the or- 
der 0.3 3 = 0.027 ~ 3%. That counting suggests that a good approximation to the NNLO 
amplitude should have an accuracy not worse than 3% which would correspond to the 
higher order not considered in the approximation. This criterion will be applied for the 
Fk/Fk quantity in the next subsection. 



Approximation II 

The above approximation is motivated by large N c behavior but some of the terms not 
included in the approximate functional (5) have a known analytical form. In practice, 
one may include those terms, i.e., the products of one-loops amplitudes from diagrams 
a and c in Fig. 1, to improve the accuracy of the approximation for certain observables. 
This second approach (called Approximation II) should also satisfy the accuracy criteria. 

Application to lattice data for F K /F n 

We now apply the analytic approximations Approximation I and Approximation 
II to the ratio Fk/Fk of meson decay constants. Fk/F k is a suited quantity for our 
exploratory exercise for several reasons: at /i = O.WGeV the genuine two-loop accounts 
less than 1% and it is known numerically ([11, 5]) allowing a direct evaluation of the 
approximation; the BMW Collaboration [12] provides 13 lattice points (considering only 
m n < 450MeV) for that ratio. 

The full NNLO calculation at physical meson masses tells us that Fx/F K — 1 = 0.14 
at @{p A ). At 0(p 6 ) there are three different contributions: tree level, 1-loopxL/, 2-loops 
and contribute to F K /F K - 1 with 0.008, 0.051 and 0.002 respectively. 

In our approximation, the 0(p 4 ) and the tree level and 1-loopxL,- at &(p 6 ) are exact 
by construction. Our estimation for the genuine 2-loops piece amounts, using M = Mk, 
—0.030 using Approximation I and —0.011 using the Approximation EL In both cases, 
the criteria of applicability is fulfilled. 

We can now proceed on fitting the 13 data points with our approximations for Fk/F k . 
We obtain for the LECs at /i = 0.77 GeV (the fitting function can be found in Ref.[l]) 
the results shown in Table 1 . 



The three parameters are strongly correlated [1]. Taking that into account, Fk/F k for 
physical meson masses is found to be 

F K /F K \ Approxl = 1.198 ±0.005, (6) 
F K /F K \ Approx n = 1.200 ±0.005, 

in good agreement with the result Fk/F k = 1.192(7)^(6)^ of Ref. [12]. Our errors 
take only into account the statistical error of the lattice values for Fk/F k but not the 
systematical once since we do not pretend to improve on the accuracy of the result but 
show the feasibility of the lattice data to explore the <ff{p b ) chiral lagrangian. 

We have include the NLO LEC L 5 in our fit because it is the only LEC that appears at 
this order and also because at NNLO, Lg is the leading piece in \/N c counting. We use 
the common fitlO [13] values for the other L ; appearing at NNLO. 

The residual dependence on the scale M is evaluated by setting M — Mk (the lattice 
value) and varying this scale by ±20%. In doing this, both F^jF n and L5 remain again 
unchanged while the LECs of <ff{p b ) vary within two standard deviations. 

Other examples 

The approximations proposed here have also been applied 2 to two more quantities. 

The first one is the ratio F k /Fq in chiral SU (3) to obtain a value for the LEC L4 and 
for the chiral parameter Fq, both badly known. The comparison of our approximations 
in that case with the numerical available results, [5] are encouraging. 

The second example is the Kft vector form factor f+(t). The particular case /+(0) 
is of interest to obtain a prediction for the CKM parameter V us . The approximations 
discussed here do not appear very promising in this case: the chiral expansion shows 
an unusual behavior and the approximations do not match with the available numerical 
results of Refs. [5, 14, 15]. 

CONCLUSIONS 

We have proposed analytic approximations for chiral SU (3) amplitudes starting from the 
structure of the generating functional of Green functions to <^( 6 ) in such a way that only 
tree-level and one-loop diagrams are needed. Our proposal provides useful formulas 
for the extrapolation of lattice data and also allows, due to their renormalization scale 
independence, for the determination of the LECs at NNLO otherwise difficult to extract 
from experimental data. These approximations are superior to the NLO amplitudes and 
to the well-known double-log approximations [16]. Concerning the l/N c expansion, 
these amplitudes contain all the leading and next-to-leading terms and all the chiral 
logarithms. This approach is specially useful in cases where the genuine two-loop 
contributions are small, in accordance to the large-Af c counting. The particular case of 
the quantity Fk/F k has this property. We fit the approximated expression of Fk/F k to 
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recent lattice data where we obtain a value for Fk/F k in agreement with the analysis 
of Ref. [12]. We also obtain a value for L 5 and for the LECs of @{p 6 ), C[ 4 + C[ 5 and 
C[ 5 + 2C[ 7 . While both Fk/F k and L5 are insensitive to the approximation made, the 
LECs of 0(p 6 ) are consistent with expectations but subject to uncertainties exceeding 
the lattice errors. All in all that suggests the potentiality to explore higher-orders on 
the chiral expansion using lattice data with an approximation that provides with more 
insight than for example a simple polynomial. 

We also suggest an Approximation II, as a modification of Approximation I, that 
includes the reducible diagrams a and c (Fig 1), when the strict Large-/V c is ignored. 
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